J/ip Polarization in 800-GeV p— Cu Interactions 
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We present measurements of the polarization of the J/ip produced in 800-GeV proton interactions 
with a copper target. Polarization of the J/ip is sensitive to the cc production and hadronization 
processes. A longitudinal polarization is observed at large xf, while at small xf the state is produced 
essentially unpolarized or slightly transversely polarized. No significant variation of the polarization 
is observed versus pr- 

PACS numbers: 13.88.+e, 14.40.Nd 



A detailed understanding of the production mecha- 
nism of charmonium is important for ongoing research at 
the Relativistic Heavy Ion Collider, where this process 
can play an important role in the search for quark-gluon 
plasma formation 0, as well as for the investigation of 
the gluon contribution to the proton spin structure 
In the non-relativistic QCD formalism (NRQCD), char- 
monium polarization can probe |4| details of the produc- 
tion process [j| that are perturbatively calculable. 

NRQCD is an effective field theory that approximates 
the full QCD Lagrangian for large quark masses. Un- 
calculable matrix elements for the production of various 
Fock-space components of the quarkonium wavefunction 
are ranked, using simple scaling rules, according to their 
order in v, the relative quark- antiquark velocity in the 
quarkonium rest frame; for charmonium, this has a value 
of ~ 0.5. The leading matrix elements can then be ex- 
tracted from fits to experimental data and used in calcu- 
lating other processes. 

A determination |(| of the various matrix elements 
from available high-energy data defines to a large ex- 
tent the expected production properties of charmonium, 
one of which is the polarization. For a spin-1 parti- 
cle, given the fraction £ of particles produced in the 
j z = ("longitudinal") state, we can define the polar- 
ization A = (1 — 3£)/(l + £); A is positive (negative) for 
transverse (longitudinal) polarization. While several of 
the intermediate cc states are color-octet states and must 
be followed by multiple gluon emission before a physical 



charmonium state is produced, heavy-quark symmetry 
implies that gluon radiation leaves the quark spins un- 
changed, providing definite predictions for the spin state 
of the final charmonium [j| . 

The polarization of the J/ip has been measured with 
relatively high statistical precision only in fixed-target 
experiments, in pion Q and proton [^, interactions 
with solid nuclear targets. No significant polarization has 
been seen in either, with the exception of an intriguing 
large longitudinal polarization at the highest xf value of 
the pion-induced data. Measurements in collider exper- 
iments, at energies and transverse momenta where the- 
oretical calculations should be more robust and which 
are free of potential complications from nuclear effects, 
suffer from low statistics [lfj|. Finally, in a recent study 
of fixed-target bottomonium production the T(15) 
was found to be largely unpolarized, while the T(25, 35) 
states had strong transverse polarization. 

In this experiment, we accumulated |l2^ a much larger 
sample — approximately 9 million — of reconstructed 
J/ip than any previous study, in interactions of an 800- 
GeV/c proton beam with a copper target. J/yb decays 
were measured using the Fermilab Meson-East dimuon 
spectrometer |l3j . which consisted of three dipole mag- 
nets, SM0, SM12, and SM3, and three stations of drift 
chambers and trigger scintillator hodoscopes. Data were 
collected during a month-long dedicated run in which 
the copper beam dump inside SM12 was used in place 
of a target; the first dipole was switched off during this 
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run. A copper absorber filtered out all hadrons and elec- 
trons produced in the interactions, allowing only muons 
to enter the spectrometer. Additional hadron rejec- 
tion was obtained using a muon identifier, consisting of 
proportional tubes and scintillator hodoscopes behind a 
thick absorber at the downstream end of the apparatus. 
Events were recorded when the trigger condition of two 
oppositely-charged muons was satisfied. 

Events were reconstructed offline from the recorded 
hits in the drift chambers. Tracks were traced back 
through the magnetic fields to the dump/target, where 
a vertex was formed, consistent with the beam position. 
Energy losses and multiple Coulomb scattering in the 
absorber and the dump were taken into account in the 
traceback. The momentum of the muons was determined 
by their bending in SM3 and this was used to calculate 
their trajectories in SM12. 

The invariant mass of the muon pair was calculated 
from the muon momenta and opening angle at the vertex. 
Since the energy loss and multiple scattering are only 
known on average and not on an event-by-event basis, 
the mass resolution, typically w 500 MeV (FWHM), was 
not sufficient to separate the J / ip and ip' peaks in this run 
using this extended target. While the two charmonium 
states cannot be separated, it is estimated that the ip' 
contributes only about 1% to the total event count, based 
on the relative production cross sections and branching 
ratios into muons. Therefore, the ip' is not considered in 
the following discussion. 

The mass distributions were plotted in bins in the 
Feynman-x variable Xf, transverse momentum px, and 
the dimuon polar angle i? in the dimuon rest frame. We 
use the Collins-Soper frame 0]. This is identical to the 
Gottfried- Jackson frame [l^, used in several earlier ex- 
periments, for pt — 0, and to a very good approximation 
equivalent even at the highest px values in this experi- 
ment. 

The invariant mass distribution in each [x f , pr)-bin 
were fitted to a Gaussian peak plus an exponential or 
polynomial background. All of the parameters, including 
the mass, were determined by the fit. The number of 
events under the peak gave the combined J/ip and ip' 
triple-differential, unnormalized cross section in Xp , Pt, 
and Distributions as a function of $ were then formed 
in each xf and px bin. 

The spectrometer and trigger acceptance was calcu- 
lated with the help of a Monte Carlo simulation of the 
J/ip production process, which included all the measured 
magnetic fields and detector efficiencies and geometry. 
The known properties of J/ip production from previous 
experiments were used. However, because of the bin- 
ning in xf and px, exact knowledge of the form of the 
cross section as a function of these two variables was 
not crucial. Events were generated with a flat d distri- 
bution, corresponding to unpolarized production. Simu- 
lated data were passed through the same analysis chain 



as the real data. The Monte Carlo reproduced quite ac- 
curately the main features of the data, including the mass 
and vertex resolutions. 

The azimuthal-angle {(p) dependence of the produc- 
tion cross section was assumed to be flat in the Monte 
Carlo, since the two interacting hadrons were unpolar- 
ized. The corresponding measured acceptance-corrected 
distribution was essentially flat, consistent with a 2% un- 
certainty in the direction of the incoming beam. The 
decay distribution was also assumed to be flat in (p, con- 
sistent with previous results Q. No attempt was made 
to extract the <j> dependence of the decay in this analysis. 
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FIG. 1: J/\& polarization parameter A versus xf in px bins. 
Solid dots are the results obtained with the 2800 A magnet 
setting, open triangles with the 2040 A setting. Only statis- 
tical errors are shown. 

The i3 distributions of the data in all bins were divided 
by the corresponding ones of the simulation, resulting in 
acceptance-corrected $ distributions. These were then 
used to calculate the polarization in xf and px bins, 
according to the formula da/dcosfl = A(l + Acos 2 $), 
where the normalization constants A were left free. Re- 
sults were obtained separately for two experimental runs 
with the current in the SM12 magnet set to 2040 and 
2800 A respectively, resulting in substantially different 
acceptances. 

Figure^shows the polarization parameter A as a func- 
tion of xf in four px bins; it can be seen that results from 
the two magnet settings are in reasonable agreement, giv- 
ing some confidence that the acceptance is understood 
and providing an estimate for the magnitude of the rel- 
evant systematic uncertainties. The xf dependence of A 



appears to be independent of px- 
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FIG. 2: J/ty polarization parameter A versus for two xf 
ranges: if < 0.45 (solid circles) and xf > 0.45 (open tri- 
angles). Statistical errors are smaller than the data points. 
Systematic errors for the small xf data are shown as a dark 
band; those for large Xf (not shown) are slightly smaller. 

Systematic errors from various sources were consid- 
ered. Inexact knowledge of the px dependence of the 
production cross section, coupled with a strong pt de- 
pendence of the acceptance versus decay angle ■&, led to 
an uncertainty of ±0.06 in A, independent of xf- Ad- 
ditional contributions included mass-peak fitting errors 
(0.04-0.08, depending on the xf bin) and uncertainties 
in the exact position (0.02) and angle (0.02-0.04) of the 
incoming beam and in the fields in the analyzing magnets 
(0.01). Uncertainties from various sources are largely un- 
corrected; they were added in quadrature for the overall 
systematic error. 

Figure |21 presents the polarization parameter A as a 
function of px for two Xf ranges, where the two data sets 
have been statistically combined. The two intervals in 
xf approximately correspond to regions were the gluon- 
gluon and quark-antiquark processes are dominant. No 
significant px dependence is seen in either region after 
the systematic errors are taken into account. At large px, 
charmonium production is understood ^(| to be domi- 
nated by gluon bremsstrahlung with subsequent fragmen- 
tation into a cc pair. In this case, the charmonium state is 
expected 0, to retain to a large degree the transverse 
polarization of the high-p^, on-shell gluon. However, the 
Pt range of the experiment, extending to 4 GeV, may not 
be sufficient to see clearly such an effect. It must be noted 
that this effect is also not observed at the high px val- 
ues available at collider energies p~£j], where, if anything, 
polarization appears to be longitudinal. In the following, 
we assume there is no significant px dependence and the 
results are presented integrated over all px- 

Figure |21 shows the polarization A as a function of 
Xf for the combined data set. The one-sigma system- 
atic uncertainty is shown as a dark band. Also shown 
are the results previously obtained by the Chicago-Iowa- 
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FIG. 3: J/ip polarization parameter A versus xf from this 
experiment (solid circles) . Statistical errors are shown as error 
bars, systematic as a dark band. Data values are tabulated 
in [Tal . The CIP @] results are also shown (open triangles). 

Princeton collaboration pj, which used a 252-GeV pion 
beam and a tungsten target. Our results suggest that the 
J I ip is produced slightly transversely polarized at small- 
to-intermediate xf- Within the systematic uncertainties 
the results are consistent with the CIP experiment pfl, 
which saw no polarization in this range, albeit with much 
larger statistical uncertainties. It should be noted that 
in this range the qq annihilation process is expected to 
play a more important role with a pion beam than with a 
proton beam, where gg fusion dominates. For xf > 0.6, 
the polarization turns to longitudinal, outside the mar- 
gins of the systematic error. This is similar to a pattern 
seen in the CIP experiment, although the behavior ap- 
pears smoother as a function of xf and begins at smaller 
xf than the sudden turn-over near xf ~ 0.85 in the lat- 
ter. However, the results are not incompatible. It is also 
interesting to note that the clear change in polarization 
in this experiment roughly coincides with the transition 
from gluon-fusion dominance to quark-annihilation dom- 
inance. 

Integrated over the entire xf range, the measured po- 
larization is A = 0.069 ± 0.004 ± 0.08 (statistical and sys- 
tematic errors). This small value, consistent with no po- 
larization if systematic errors are taken into accountjis in 
agreement with previous proton-beam experiments 0,0, 
which had insufficient statistics to study the xf depen- 
dence. 

As a cross-check of the analysis, the polarization of 
the Drell-Yan continuum was also studied using the same 
technique. This can be done only for dimuon invariant 
masses greater than 4 GeV, since the J/ip peak domi- 
nates the spectrum at lower masses. An additional com- 
plication arises from random coincidences of uncorrelated 
pairs of muons, which in the case of the J /ip were re- 
moved by the fitting procedure. These were estimated 
and subtracted from the Drell-Yan sample by studying 
the distributions of same-charge muon pairs. A fit to 
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the formula 1 + A cos 2 d for dimuon masses from 4 to 7 
GeV gave A = 0.98 ± 0.04, in good agreement with the 
100% transverse polarization expected and previously ob- 
served ^lj m the Drell-Yan process. While the kinematic 
range is not the same as for the J/ip measurement, this 
agreement with expectations increases confidence in the 
soundness of the J/ip analysis. The level of statistics 
were not adequate to investigate any xf dependence of 
the Drell-Yan polarization. 

The results are not in agreement with published pre- 
dictions for J/ip polarization based on NRQCD, which 
are in the r ang e 0.31 < A < 0.63 or with similar 
predictions [19j | based on an early model that considers 
only color-singlet intermediate states. The small positive 
values at xf ^ 0.5 are very similar to the results obtained 
for T(IS') by this experiment in a similar px range, and in 
sharp contrast to the essentially 100% transverse polar- 
ization of the (unresolved) T(25) and T(35*) states [Tl| . 

To understand the behavior of the polarization, all 
sources of J/ip production must be considered. While 
6-quark production is not a major source at these en- 
ergies, almost half of the produced J/ip's are the decay 
products of higher-mass charmonium resonances, mainly 
the 3 Pj states Xcj- Feed-down from \c2, produced exclu- 
sively in the J z = 2 state from gluon-gluon fusion, results 
in 100% transverse J/i[>'s, increasing the observed values 
of A. The turn-over at high xf may reflect the transi- 
tion, at xf — 0.6, to the quark-annihilation graph |20| . 
which produces a mix of Xc2 spin states. As for the overall 
level, it cannot be explained without a substantial contri- 
bution from XcXi which in gene ral p roduces longitudinal 
J/ip's. Beneke and Rothstein [18J originally estimated 
the Xci contribution to be about 10 times smaller than 
that from Xc2 , but it was later found experimentally |2l| 
to be of comparable size. NRQCD calculations which 
retain higher orders in v can also accommodate substan- 
tial Xci contributions [22l |. It would probably require a 
contribution near the experimental and theoretical upper 
limits in order to obtain A values as small as experimen- 
tally measured, within the context of NRQCD. 

In addition, it must be noted that no existing calcu- 
lation takes into account nuclear effects, which strongly 
affect the production cross section [2j|. The formation 
length of the J/ip at these energies is generally longer 
than the nuclear size jiij and it is conceivable that color- 
singlet and -octet components of the wavefunction are 
absorbed differently while propagating through the nu- 
clear medium, resulting in a different mix of Fock states 
compared to the free-nucleon production, thus altering 
the polarization. Furthermore, the polarization of the Xc 
states themselves is predicted to have a nuclear depen- 
dence which would feed down to the J/ip polariza- 
tion. Finally, at the highest values of xf, higher- twist 



effects may become important [l9T | . 
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